Spontaneous microcavity-polariton coherence across the parametric threshold: 

Quantum Monte Carlo studies 
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We investigate the appearance of spontaneous coherence in the parametric emission from planar 
semiconductor microcavities in the strong coupling regime. Calculations are performed by means 
of a Quantum Monte Carlo technique based on the Wigner representation of the coupled exciton 
and cavity-photon fields. The numerical results are interpreted in terms of a non-equilibrium phase 
transition occurring at the parametric oscillation threshold: below the threshold, the signal emission 
is incoherent, and both the first and the second-order coherence functions have a finite correlation 
length which becomes macroscopic as the threshold is approached. Above the threshold, the emission 
is instead phase-coherent over the whole two-dimensional sample and intensity fluctuations are 
suppressed. Similar calculations for quasi-one-dimensional microcavities show that in this case the 
phase-coherence of the signal emission has a finite extension even above the threshold, while intensity 
fluctuations are suppressed. 

PACS numbers: 71.36.+C, 42.50.Ar, 42.65.Yj, 02.70.Uu 



I. INTRODUCTION 

The recent progresses in the growth and manipulation 
techniques of semiconductor heterostructures have re- 
cently led to the realization of planar semiconductor mi- 
crocavities in which a cavity-photonic mode is strongly 
coupled to an excitonic transition, so that the resulting 
elementary excitations consist of linear superpositions of 
cavity-photon and exciton, the so-called polaritonsi^^. 

The polaritons satisfy the Bose statistics, at least in a 
low-density regime, and have recently started to be con- 
sidered as potential candidates for the study of many- 
body effects, e.g. superfluidity-, in regimes completely 
different from the ones usually considered in experiments 
with Helium fluids 5 or ultracold atomic gases 6 . 

A further advantage of polaritonic systems consists of the 
fact that polaritons can be optically injected in the cav- 
ity by simply shining an incident light beam onto it, and 
their properties can then be inferred from the observation 
of the emitted light. In particular, the coherence prop- 
erties of the light emission closely reproduce the ones of 
the in-cavity polaritonic field!. 

Thanks to their excitonic component, polaritons have 
strong binary interactions, which have been exploited 
to obtain spectacular polariton amplification and opti- 
cal parametric oscillation effects^&iLiil. Differently from 
conventional optical resonators with a discrete set of op- 
tical modes, planar microcavities are spatially extended 
systems with a continuum of transverse modes and the 
parametric oscillation process automatically selects the 
most favourable one in which to generate the coherent 
signal field. 

Simultaneously with these experiments, an active debate 
has taken place in the community about the relation be- 



tween the appearance of spontaneous coherence in op- 
tical parametric oscillation in planar microcavities and 
Bose-Einstein condensation of excitonsii. The present 
paper is devoted to shine light on analogies and differ- 
ences between the two phenomena: both of them are in 
fact characterized by the appearance of coherence and 
long-range order in a Bose HeldS^^. There is however 
a major difference: while Bose-Einstein condensation is 
an equilibrium phenomenon, parametric oscillation takes 
place in intrinsically out-of-equilibrium systems, where 
the stationary state arises from a dynamical equilibrium 
between the driving by the pump and the dissipative ef- 
fects due to losses. From this point of view, the para- 
metric oscillation threshold is closely related to the laser 
thresholdi^i 5 ", although standard laser cavities have a fi- 
nite number of discrete modes. 

In order to perform quantitative calculations, a Monte 
Carlo technique based on the Wigner representation of 
a quantum fieldi^ has been generalized to the case of a 
pair of coupled excitonic and cavity-photonic fields and 
the resulting numerical code has been used to calculate 
the one-time correlation functions of the in-cavity field in 
the stationary state, in particular the first- and second- 
order ones. As already mentioned, these univocally de- 
termine the spatial coherence properties of the emitted 
light which can be measured in the experiments. 

The advantages of the Wigner-Monte Carlo approach 
with respect to previous works on the subject are many- 
fold: first of all, we keep track in a complete way of 
the fluctuations of the polaritonic field, even in the crit- 
ical region around the threshold, where fluctuations are 
large and the linearized approaches around the mean- 
g e i ( j9 i i8 4 i9 fgji Furthermore, no few-mode approxima- 
tioni2i22i2ii2£ is here performed: the parametric emission 
into all available modes is automatically taken into ac- 
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count, as well as the correlations between all modes. 

Above threshold, the signal emission in planar microcav- 
ity geometries occurs in the single, most favourable, mode 
(which the numerical code is able to automatically select) 
and therefore shows the typical spatial coherence proper- 
ties of a Bose-Einstein condensate. As it is not inherited 
from the pump beam, this coherence is a spontaneous 
one. 

On the other hand, the parametric emission from 
reduced-dimensionality systems such as photon wires 2 ^ 
shows dramatically different properties: as long- 
wavelength fluctuations of the signal phase are here able 
to destroy the long-range spatial coherence, the thresh- 
old is no longer well defined and is replaced by a smooth 
crossover where a local order sets in and the intensity 
fluctuations of the signal emission are suppressed. Such 
a behaviour is closely related to the so-called quasi- 
condensation phenomenon of equilibrium Bose gases in 
reduced dimensionalities 2 ^ 2 ^*^. 

Related questions about the spontaneous appearance of 
long-range coherence and the effect of a reduced dimen- 
sionality have been recently adressed2L2& for the simpler 
case of a parametric oscillator formed by a planar res- 
onator containing a slab of generic x^ medium without 
any excitonic resonance. 

The plan of the present paper is the following: in seclITI 
we introduce the field-theoretical Hamiltonian of the sys- 
tem and we relate the properties of the emitted light 
to the ones of the in-cavity polaritonic field. In sec lIIII 
we work out the corresponding mean-field theory and 
we use it to identify the parameter range which is most 
suited for the following analysis. In sec II VI we review the 
Wigner representation formalism and we generalize it to 
the case of coupled cavity-photon and exciton fields. We 
also describe the principles of the numerical code used to 
obtain the results which are discussed in the following. 
In sec0 the central numerical results of the paper are 
shown: the behaviour of the first- and second-order co- 
herence properties of the signal emission are presented for 
pump parameters spanning across the parametric thresh- 
old. These results are the basis of the discussion of the 
following sec I VII where analogies and differences between 
our driven-dissipative system and the Bose-Einstein con- 
densation phenomenon of equilibrium statistical mechan- 
ics are pointed out. Sec lVllI an alternative interpretation 
of the parametric threshold in terms of modulation in- 
stability and spontaneous breaking of the translational 
symmetry is given and its consequences on the intensity 
correlation function are put forward. The effect of a re- 
duced dimensionality is discussed in sec lVIIII where the 
link with the concept of quasi-condensate is evidentiated. 
Finally, in sec lIXI we give our conclusions, and we sketch 
the perspectives of the work. 



II. THE SYSTEM HAMILTONIAN 

A model which is commonly used to describe a pla- 
nar microcavity with a quantum well excitonic resonance 
strongly coupled to the cavity mode is based on the field- 
theoretical Hamiltonianii: 



n 



Jdx J2 *i « K- + Vi (x) 5ij ] * j (x) 



zj={X,C} 



dxh/3i n cE p (x,t)¥ c (x) + h.c. , (1) 



where x is the in-plane spatial position. The field opera- 
tors \&x,c( x ) respectively describe excitons (X) and cav- 
ity photons (C) and, provided the exciton density is much 



lower than the saturation density n s 
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they satisfy 



Bose commutation rules, [^(x), $t(x')] = <5 2 (x — x') 6ij. 
For the sake of simplicity, we consider here a single polar- 
ization state, but a generalization of to include sev- 
eral polarization states and describe the spin dynamics 
discussed e.g. in 

ReflU is straightforward. 
The linear Hamiltonian h° is: 



h u = h 



wx(-*V) 



wo(-*V) 



(2) 



where wc(k) = ui^ + k 2 /fc| is the cavity dispersion 
as a function of the in-plane wavevector k and k z is the 
quantized photon wavevector in the growth direction. 
The exciton dispersion is instead much weaker and can 
be well approximated by a constant w^(k) = lux- The 
quantity flu is the Rabi frequency of the exciton-cavity 
photon coupling: the eigenmodes of the linear Hamilto- 
nian © are therefore linear combination of the cavity 
photonic and excitonic modes. They are called upper 
(UP) and lower (LP) polariton and their dispersion is 
given by u L p(jup) (k). 

The nonlinear interaction term is due to exciton-exciton 
collisional interactions and can be modelled by a repul- 
sive (g > 0) contact potential. The anharmonic exciton- 
photon coupling has a negligible effect in the regime con- 
sidered in the present study-i 7 -. Vx,c( x ) are external po- 
tential terms acting on the excitonic and photonic fields 
and can be used to model etched cavities 23 . 

Differently from other Bose systems such as ultracold 
trapped atoms or liquid Helium, our polaritonic system 
is intrinsically out of equilibrium. Polaritons are in fact 
injected in the cavity by incident light beams and decay 
through several dissipative channels. In the present pa- 
per, polaritons are coherently injected in the cavity by a 
coherent pump laser, which is taken into account in the 
Hamiltonian Q by the driving term proportional to the 
incident electric field E p (x,t) of the pump. For a given 
cavity, the coefficient (3i nc can be calculated from the re- 
flectivity of the front mirror through which the cavity 
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is excited^. Both the spatial and the temporal depen- 
dences of E p (x,t) can be experimentally controlled, so 
to obtain a variety of different temporal (e.g. monochro- 
matic, bichromatic, pulsed, etc.) and spatial (e.g. plane 
wave, gaussian, optical vortex, etc.) patterns^. 

Dissipative effects lead to a damping of both the excitonic 
and cavity photonic fields. Using the standard quantum 
theory of damping, they can be modelled by additional 
terms in the master equation giving the time-evolution 
of the density matrix p of the systeml&: 



dp 
~dt 



■[n, P ] 



dk 

2^ 



£ 7i (k)[^(k)p*J(k) 



={X,C} 



--($t(k) h(k-) P + P*m *i(k)) 



(3) 



The operators $x c(k) are the Fourier transform of the 
field operators i&x,c(' x )- In the present paper, we shall 
consider the case of momentum-independent jc.x which 
summarize dissipative effects of different origin, e.g. exci- 
tonic scattering and non-radiative recombination, back- 
ground absorption by the material forming the structure 
as well as radiative emission into the continuum of modes 
outside the cavity. 

From the experimental point of view, information on the 
state of the polariton field in the cavity is generally re- 
trieved by detecting and characterizing the light emission 
from the cavity through either the front or the back mir- 
ror. In particular, the normally-ordered expectation val- 
ues for the light emitted through the back mirror (called 
in the following transmitted light )a* are proportional to 
the corresponding ones for the photonic component of 
the in-cavity polaritonic field. For instance, the near- 
field amplitude of the transmitted field is 

(^(x)£ tr (x')) = |/3, r | 2 (^(x)* c (x')), (4 ) 

where the coefficient j3t r relating the transmitted field to 
the in-cavity one has to be calculated from the transmit- 
tivity of the back mirror—. Similar formulas hold in the 
far-field. The far-field angular pattern is proportional to 
the in-cavity occupation number: 

(Ei(k)E tr (k)) = \f3 tr \ 2 n c (k) = 

H/3 tr | 2 (^(k)^ c (k)) (5) 

and analogous expressions hold for its higher-order, 
normally-ordered moments such as: 

(El(k) Ej r (k') E tr (k') E tr (k)) = 

= |/3 tr | 4 (*t,(k) %(k>) $ c (k') $ c (k)>. (6) 

This procedure has recently been used to experimentally 
determine the first- and second-order correlation func- 
tions of the polaritonic field by measuring the coherence 
and noise properties of the emitted light 7 . 



III. MEAN-FIELD THEORY 

A simple way of coping with the field Hamiltonian 
is to study the corresponding classical field theory once 
the quantum field operators \I/x,c( x ) are replaced with 
classical C-number fields -0x,c( x ) defined as the expec- 
tation value of the quantum field ?/>x,c( x ) = (^x,c( x ) )■ 
For a system of weakly-interacting bosons like microcav- 
ity polaritons, fluctuations are generally weak provided 
the mean-field solution is dynamically stable, so that the 
mean-field solution provides accurate predictions for the 
coherent part of the polariton field, and therefore the 
coherent emission from the microcavity. This includes 
e.g. resonant Rayleigh scattering in the presence of de- 
fects^. In the present paper, the mean-field prediction 
for the threshold intensity will be used as a guideline 
in the choice of parameters to be used in the numerical 
calculations. 

As we are most interested in the basic properties of the 
parametric instability threshold, we have tried to iden- 
tify a configuration where the intrinsic phenomenology 
is not masked by additional effects. In particular, we 
have considered the case of a homogeneous cavity with 
no external potential Vx,c = and a transversally homo- 
geneous, monochromatic and cw laser field of frequency 
uj p which incides on the cavity at an angle p with the 
normal to the cavity plane. Provided the beam actually 
used in the experiment is wide enough in the transverse 
direction, this approximation is expected to be accurate. 
The pump electric field has then the simple form: 



E p (x,t) =E p e i ( k >- x -"* t) . 



(7) 



the wavevector k p being fixed by the incidence angle 
through k p — sm8 p lo p /c. The incident intensity I p is 
proportional to |-E P | 2 : 



I P = ^ \E P \ 2 . 



(8) 



Given the translational symmetry of the system, we can 
look for stationary solutions in the same plane- wave form 
ip x ,c(x,t) = V>f,c e i( - k " x -^ as the incident field. The 
coefficients ipx c can De obtained from the equation of 
state: 



( Wjf (kp) ~ " P - l - lx + 9 \n I 2 ) Vx + Mr rc - (9) 
(coc(kp) -u P - l --ic) Vc + ^Ri>x = -PincE p (10) 



and the stability of the solution has to be verified by 
looking at the eigenvalues of the Bogoliubov matrix 4 : 
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/ uj x + 2g \j) 



C 



x 

n, 



»7x 
2 



ss * 2 „ — 2'ik D x 








2 



5 ^2 e 2ik p x o ^ 



uj x +2g\r x s \ 2 )+u p - l ^f- -0 R 

-w c (-iV)+w p -f / 
(11) 



If the imaginary parts of all the eigenvalues ui^p LP {\&) 
of the Bogoliubov matrix C are negative for all k, the 
mean-field solution is dynamically stable; if at least one 
of them is positive, a different, spatially structured, sta- 
tionary solution has to be found, e.g. with a sinusoidal 
intensity profile, and its stability has again to be veri- 
fied by means of the appropriate Bogoliubov matrix. An 
analytical study of the mean-field stationary state in the 
parametric oscillation regime when translational symme- 
try is broken is beyond the scope of the present work. 
Simple discussions under a simple three-mode approxi- 
mation taking into account only the pump, the signal 
and the idler modes can be found e.g. in Ref ll3l20ll2lll2S 
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FIG. 1: Real (left) and imaginary (right) parts of the Bogoli- 
ubov dispersion for different values of the pump frequency 
hS p = h(u p -WLp(k. P )) = -0.65 meV (a,b), -0.54meV (c,d), 
— 0.49meV (e,f). At the mean-field threshold hS p = fi<5p C ' « 
— 0.52meV, the mean-field energy hg\tpx\ 2 — 0.13 meV. For 
the sake of clarity, only the branches corresponding to the 
LP have been traced, the UP ones being far away in energy. 
Cavity parameters: Hu>q — hujx = 1.4 eV, 2MIr = 5meV, 
h-yx = 0.15, meV, fvyc = 0.24 meV. The lower polariton en- 
ergy at linear regime is huuLp(k. P ) = 1.39935 eV. 

In order to keep the physics the simplest, we have cho- 
sen a fixed value of the pump intensity I p and we have 
spanned its frequency uj p from far below the lower po- 
lariton frequency at linear regime WLp(k p ) in the up- 
ward direction. The incidence angle 9 p has been taken 
slightly larger than the so-called magic angle, so that 
2uj L p(k p ) > uj l (0) + uj LP (2k p ). The corresponding Bo- 
goliubov dispersions are plotted in fig^ as expected, 
the first feature that takes place is the parametric insta- 



bility when the ± branches, which are coupled by the 
anti-hermitian terms of l|llf) . touch around k«0. This 
condition 2uj p — uj^ p (0) + w^p(2k p ) can be interpreted 
as the parametric process (k p ,k p ) — > (0, 2k p ) becoming 
resonant on the final state. 

For our specific choice of 9 p and u> p , notice that this hap- 

(c) 

pens at a value to p — u> p still red-detuned with respect 
to LULp(h p ). This guarantees that we arc still in an opti- 
cal limiter regime and no single-mode optical bistability 
effects— has yet taken place. In the following of the paper, 
the frequency u> p c ^ will be designed as the parametric in- 
stability threshold and the wording below/above thresh- 
old will refer to u> p respectively red- or blue-detuned with 

(c) 

respect to u> p . This definition of the threshold in terms 
of frequency is alternative but equivalent to the standard 
one in terms of the intensity and allows us to avoid those 
complicate additional effects that have been pointed out 
in Refl22I For notational simplicity, we shall also intro- 
duce the quantity 6 P = lo p — WLp(k p ); at the threshold 

s p = 4 c) < o. 



IV. THE WIGNER MONTE CARLO METHOD 

A more complete analysis of the parametric threshold can 
be carried out by numerically solving the stochastic equa- 
tions of the Wigner representation of the quantum Bosc 
field. This technique allows one to study the dynamics of 
the field taking fully into account the multi-mode struc- 
ture of the field and eventually the spontaneous breaking 
of the translational symmetry, as well as the fluctuations 
around the mean-field solution, however large they hap- 
pen to be. In the present section, we shall review the 
formalism of the Wigner representation of the quantum 
Bose field as applied to a pair of coupled photonic and 
excitonic fields with a quartic interaction term. Recently, 
a similar Wigner approach has been applied to study the 
dynamical properties of ultracold atomic gases^i. 

The Wigner representation of the two-component quan- 
tum Bose field ^(x) (i = {X, C}) consists in a quasi- 
probability distribution function W[^/)j(x)] in the func- 
tional space of the C-number fields ipi (x)i&. Any numer- 
ical implementation requires one to work on a discrete 
and finite spatial grid: in our simulations, a rectangular 
grid of Af x x Ny points with uniform grid spacings £ x<y 
along respectively the x, y directions has been used, with 
periodic boundary conditions. The integration box has 



5 



therefore sides L x>y = Af x ,y^x,y 

Under a diluteness condition that can be written as 
7i ^ 9/dV (dV — £ x £ y being the cell volume), the 
evolution of the quasi-probability distribution W can 
be approximated in the so-called truncated Wigner ap- 
proximatior^^ by a true Fokker-Planck equation with 



a positive-definite diffusion matrix and without higher- 
order derivative terms, which guarantees the positivity of 
W at all times and the possibility of mapping the prob- 
lem onto stochastic differential equations for the fields 








Fx(x)+%(hMx)| 2 -^ 



2 




^dW c {*) 



dt- 



(12) 



where dWx,c(x)> are zero-mean, independent, complex- 
valued, white noise terms such that 



cW{(x)dWj(x') = 



(13) 

dWi(x)dW;(x') = 2cft<5 x , x '%. (14) 

In the most relevant case of an initially empty microcav- 
ity, the initial fields at t = are independent zero-mean 
complex Gaussian variables such that: 



t/ji(x)ipj(x) = 



^(x)V 3 *(x) = YKv SiJ 



(15) 
(16) 



In a practical Wigner-Monte Carlo simulation, the 
Wigner distribution of the initial state at t = is sam- 
pled with a large number J\f of independent realizations 
(the data shown in the present paper have been calcu- 
lated from a number J\f of realizations of the order of a 
few thousands). They are then let evolve according to 
the stochastic equation l|12f) with a random noise which 
is obviously different for each realization. 

The expectation value of any observables at time t is 
obtained as an average over the different realizations of 
the corresponding classical quantity, paying attention to 
the fact that the moments of the Wigner function give 
the expectation value of totally symmetrized operators, 



e.g. 



1 



{O x . . . N ) W = — £ (0 P{1) . . . P{N) ), (17) 
' p 

where the sum is taken over all the permutations P of 
N objects, each of the O a 's stands for some field com- 
ponent ipi (x) or ip* (x) and O a is the corresponding field 
operator. As a simplest example, 

<V£(x) 1>ctf)) w = \ <*LW*c(x') + *cV)*U*)) 

(18) 



In the following of the paper, we shall use this Wigner- 
Monte Carlo method to numerically study the one-time 
correlation functions of the in-cavity fields in the station- 
ary state. These are obtained by letting the evolution to 
run for a time interval longer than all the characteristic 
time scales of the problem, until the stationary state is 
attained (in practice an evolution time of at most 5 ns is 
used). At this point, the expectation values of the observ- 
ables are evaluated from the corresponding averages over 
the Wigner distribution. Note that we have chosen the 
excitation parameters in such a way that the pump has a 
simple optical limiter behaviour and complex bistability 
effects are avoided. 

For instance, the first-order correlation function in the 
near-field has the form: 



<^(x)* c (x'))= (V£(x)Vc(x')> 



<5x,x' 



w 2AV 



(19) 



while the second-order one reads: 



(%x)^(x')*c(x')*c(x)> = 
|^(x)| 2 |Vc(x')| 2 ' 



2 AV 



(1 + Sx.x') 



(hMx)| 2 + KMx')l : 



1 



2AVJ/W 



(20) 



Similar formulas hold for the k-space observables which 
describe the emission pattern in the far-field. 

The study of multi-time correlation functions which are 
involved in the temporal coherence properties requires 
the calculation of non-trivial different-time commutators 
or a specific treatment of quantum noise incident on the 
cavity 2 ^. All these issues go beyond the scope of the 
present work and will be the subject of a forthcoming 
publication. 
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V. COHERENCE PROPERTIES OF THE 
PARAMETRIC EMISSION 



The Wigner-Monte Carlo method reviewed in the previ- 
ous section can be applied to study the parametric emis- 
sion for Lu p spanning across the threshold. Depending on 

(c) (c) 

whether lo p < u> p or w p > u>p , the emission shows in 
fact dramatically different features which appear clearly 
in the angular intensity pattern of the far-field emission 
as well as in its first- and second-order coherence prop- 
erties. 



A. Below the parametric instability threshold 

1. Far-field angular pattern 

Below the threshold, the far-field emission patterns 
shown in Fig|3 are characterized by a strong and nar- 
row peak corresponding to the mode at k p being di- 
rectly pumped by the incident laser field, and weaker and 
broader peaks corresponding to the luminescence from 
other modes due to parametric processes (k p ,k p ) — ► 
(k s ,k,; = 2k p — k s ). As the pump frequency lu p is in- 

(c) 

creased and approaches the threshold at w p , not only 
does the intensity of the luminescence peaks increase, 
but it is also spectrally narrowed in k space, until it be- 
comes a very narrow peak in the vicinity of the instability 
threshold. 

As discussed in Refl33l and experimentally verified in 
Refl34l the parametric luminescence has remarkable cor- 
relations properties, in particular a significant degree of 
multimode squeezing between the signal and idler. How- 
ever, as these properties do not play a central role in 
the coherence properties of the signal alone, we shall not 
analyze them in detail in the present paper. 



2. First- and second-order coherence 

The real space counterpart of the linewidth reduction is 
an increase of the coherence length. In the experiments 7 , 
a region S in k-space around the peak of the signal emis- 
sion (the region inside the squares in fig|2J) is selected 
by means of suitable combination of diaphragms. The 
corresponding real-space correlation functions are then 
defined as: 



G (1) (x,y) 
G( 2 >(x,y) 



* c (x)* c ( y ; 



(21) 



C(x)C(y)*c(y)*c(x)). (22) 



in terms of the field ^>c( x ) operators 

1 



*c(k)e 



ikx 



(23) 




o o 

O -0 

O 



-0.2 0.2 

k x (urn" 1 ) 



FIG. 2: Far-field emission pattern for increasing values of 
the pump frequency h8 p = — 0.55meV (a), — 0.53meV (b) be- 
low the mean-field threshold, and —0.52 meV (c) close to the 
mean-field threshold. The pump mode at k p by far saturates 
the grey scale, while the idler emission is hardly visible on 
this scale. The solid rectangles indicate the k-space region 
which forms the signal emission. Throughout the whole pa- 
per, the correlation functions of the selected signal emission 
are denoted with a bar (e.g. and <^ 1,2 '). Panel (d): cut 

for k y — of the far-field emission pattern shown in the pan- 
els (a-c). For clarity reasons, only the part of the spectrum 
corresponding to the signal emission is shown. Same cavity 
and pump parameters as in figQ physical size of the system 
L x = L y = 25 /im. 



where the sum over k vectors is here restricted to the 
selected region S only. The bar over G^ 1 ' 2 ) indicate that 
these correlation functions are restricted to the selected 
signal emission only. As usual, the translational invari- 
ance of the system guarantees that all two point correla- 
tion functions depend only on the distance x — y. 

In figElwe have plotted the first order coherence function 
for the signal emission: 



5 (1) (x) 



G«( 



G«(y,y) 



(24) 



For uj p well below the threshold, g^\x) goes from 1 to 
within a transverse coherence length £ c which due the 
Wiener-Khintchine theorem is proportional to the inverse 
of the k-space linewidth. As the threshold is approached, 
l c increases. At the threshold, g^(x) is everywhere fi- 
nite so that the coherence extends over the whole (finite) 
sample. From an experimental point of view, (^^(x) can 
be measured by first angularly filtering out the emission 
in the k-space region S centered around k s , and then 
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FIG. 3: First-order coherence function g^{x.) of the selected 
(see solid rectangle in fig® signal emission for values of the 
pump frequency hS p — — 0.55meV (dotted), hS p — — 0.53meV 
(dashed) and hS p — —0.52 meV (solid) respectively below and 
around the mean-field threshold. The cavity and the other 
pump parameters are the same as in fig0 



experiments like the ones reported in Ref0- As pre- 
viously, the bar over g@> indicates that this coherence 
function includes the selected signal emission only. 

For lo p well below the threshold Wp , g( 2 )(x) shows the 
typical Hanbury-Brown and Twiss bunching^ at short 
distances |x| <C t c , where its value is close to 2, and then 
it tends to 1 for |x| i c . This behaviour can be eas- 
ily understood in terms of the Gaussian structure of the 
field: well below the threshold, all modes of the coupled 
excitonic and photonic fields in the microcavity are in a 
Gaussian state, exception made for the one at k p directly 
driven by the pump. This means that Wick theorem ap- 
plies and higher-order correlation functions can be conse- 
quently factorized into products of first-order onesS*. As 
the parametric emission gives non-vanishing anomalous 
correlations only between the signal and the idler modes, 
all the anomalous correlations appearing in the Wick fac- 
torization of the signal second-order correlation function 
(x, y) vanish and this can be factorized as: 

( X , y) ~ G« (x, x) G« (y, y) + | G« (x, y) | 2 (26) 

so that g^ (x) is 



performing a sort of Young's two-slit experiment like the 
one in Refl7l 



>(x)~l+|ffW(x)| 



(27) 




x (|im) 



FIG. 4: Second-order coherence function g^(jx.) of the se- 
lected (see solid rectangle in fig® signal emission for values 
of the pump frequency h8 p = —0.53 meV (dashed) —0.52 meV 
(solid) respectively below and around the mean-field thresh- 
old. The cavity and the other pump parameters are the same 
as in fig0 



The second-order coherence function of the signal emis- 
sion g^ is defined as: 



5 (2) W 



(25) 



[G«(y, y )p 

and experimentally can be accessed by noise correlation 



This form is consistent with the numerical results shown 
in Figs® and® In particular, the coherence length of 
G' 2 ' is roughly half the one of . Note that the decom- 
position i|2b|) would not be valid for quantities without a 
bar, e.g. G^ 2 ^(x) that will be introduced in (|29[l) since in 
that case the anomalous correlation between signal and 
idler gives an important contribution. 

Closer to the threshold, the field is no longer in a Gaus- 
sian state: the short-distance value g^(0) which de- 
scribes the intensity fluctuations of the signal emission 
is reduced to a value close to 1 as it happens for a coher- 
ent field. 



B. Above the parametric instability threshold 



Far-field angular pattern 



As the pump frequency lo v goes through the critical value 

(c) 

ujp , the signal and idler emission concentrates in a sin- 
gle signal/idler pair of modes at k Sj i and their intensity 
becomes of the same order as the one in the pump mode 
at k p (Fig®. 

In Fig® we have plotted the intensity in the pump and 
the signal modes as a function of the pump frequency 
Lu p . While the pump intensity N p is simply the occupa- 
tion number ncQ&p) of the mode at k p , the signal and 
idler intensities N s i are defined as the sum of the occu- 
pation numbers rip (k) for k within the regions contained 
in respectively the solid and dashed rectangles in Fig® 
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Both the pump and the signal intensities show a thresh- 

(c) 

old behaviour at ujp : the pump intensity N p mono- 
tonically increases upto the threshold, and then starts 
slightly decreasing because of the presence of the finite 
signal [figEJa)]. The signal one N s is very weak be- 
low the threshold, and has a sudden increase close to 
Wp C ' [figEJb)]. As compared to the mean- field predic- 
tions^, the transition is rounded by quantum fluctua- 
tions, which give indeed a finite signal intensity also be- 
low the threshold and smoothen the threshold. Note that 
fluctuations are also responsible for the slight shift of the 
threshold which can be observed in the figure. The be- 
haviour of the idler intensity (not shown) is strictly anal- 
ogous to the one of the signal one. 
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FIG. 5: Upper (a) panel: far-field emission pattern for 
a pump frequency hS p — —0.5 meV (above the mean-field 
threshold). The solid (dashed) rectangles indicate the k-space 
region which forms the signal (idler) emission. Throughout 
the whole paper, the correlation functions of the selected sig- 
nal emission are denoted with a bar (e.g. G' 2 ^ and <^ 1,2 '). 
Lower (b) panel: cut for k y — of the far-field emission pat- 
tern shown in panel (a). The cavity and the other pump 
parameters are the same as in fig0 

Far above the threshold, additional, weaker, peaks ap- 
pear due to multiple scattering process^, but they are 
filtered out when the signal emission is selected. A dis- 
cussion of these coherent peaks, as well as of the novel 
spectral and angular features shown by the weak inco- 
herent luminescence when we are far above the threshold 
go beyond the scope of the present paper and will be 
postponed to future work. 



2. First- and second-order coherence function 



As in sec IV A 2l we have selected the region S in wavevec- 
tor space corresponding to the signal emission (within 
the solid rectanble in Fig|5J) and we have determined its 
coherence properties. 



2300 1 — 
a. 2250 - 
2200 - 
2150- 
2100 - 
2050 B 



0.005 0.01 0.015 0.02 



co -co 10 ' [meV] 



o- (b) 



E 




co -co 1 "' [meV] 
P P 



FIG. 6: Intensity in the pump (a) and signal (b) modes as 
a function of the pump frequency uj p . Solid lines: Wigner- 
Monte Carlo results. Dashed lines: mean-field prediction. 
The dotted line in panel (a) gives the intensity in the pump 
mode for the unstable homogeneous mean-field state of 
llOH . The cavity and the other pump parameters are the same 
as in fig'0 
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FIG. 7: First- (dashed) and second-order (solid) coher- 
ence functions g^ 1,2 \x.) of the selected (see solid rectangle in 
figHJ signal emission for a pump frequency hS^ — — 0.5meV 
(above the mean-field threshold). The cavity and the other 
pump parameters are the same as in fig0 



The first- and second-order coherence functions 5*- 1,2 ^(x) 
are plotted in figd 5^H X ) i s nearly flat and equal 
to 1, which indicates that the first-order coherence ex- 
tends over the whole two-dimensional sample, or, in other 
terms, that the emission is phase coherent. As expected 
for a coherent field, the second-order coherence function 
<7( 2 )(x) is also almost flat and its value close to 1: in phys- 
ical terms, this means that the intensity fluctuations of 
the signal emission are weak. 
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Both these qualitative behaviours are in agreement with 
the experimental observations reported in Ref 7. 



VI. ANALOGY WITH THE BOSE GAS AT 
THERMAL EQUILIBRIUM 

The analogy between the threshold of a optical paramet- 
ric oscillator or of a laser and a sort of non-equilibrium 
second-order phase transition was put forward in the 
'70 o 14 i 15 , but its study was mostly focussed on the case 
of few-mode systems, such as optical resonators with dis- 
crete modes. Its generalization to spatially extended sys- 
tems showing a continuum of modes is instead a much 
more recent topic from both the theoretical^! and the 
experimental 7 points of view. In the present section, we 
shall try to interpret the results presented in the previ- 
ous section in terms of a sort of non-equilibrium Bose- 
Einstein condensation. 

At the level of mean-field theory, the analogy is immedi- 
ate: in usual statistical mechanics, the onset of a phase 
transition is signalled by the thermodynamical instability 
of the symmetric state and the consequent appearance of 
a non- vanishing value of the order parameter in the equi- 
librium state. For the specific case of the Bose-Einstein 
condensation transition, the order parameter consists of 
the expectation value of the Bose field ^(x). 

In the present non-equilibrium case, the mean-field 
threshold corresponds to the point where the homoge- 
neous solution © and (fTUf> becomes unstable and a finite 
signal and idler amplitudes appear in the new stable sta- 
tionary state. As in the standard theory of the optical 
parametric oscillator, the absolute phases of the signal 
and the idler are not determined. The Hamiltonian is 
in fact invariant under the U(l) symmetry which sends 
ip(s,i) ^,(«,») g±j0 anc j the specific value of the phase is 
randomly chosen at each realization by the spontaneous 
symmetry breaking mechanism. The chosen value of the 
order parameter is the same throughout the whole sys- 
tem. 

Quantitative results to illustrate the physics of this 
spontaneous symmetry breaking are obtained from the 
Wigner-Monte Carlo calculations including the fluctua- 
tions around the mean-field which smoothen the thresh- 
old. Because of the fluctuations, the signal and idler 
amplitudes no longer vanish below threshold, but as the 
signal and idler emissions take place into many available 
modes which are close to resonance, they are incoherent 
both spatially and temporally. As the threshold is ap- 
proached, the characteristic coherence length which de- 
termines the spatial behaviour of the first- and second- 
order coherence functions g^ 1,2 '(x) increases and becomes 
macroscopic at the threshold. Simultaneously, the inten- 
sity fluctuations, quantified by g^(0), are quenched. 

Above the threshold, the parametric emission is concen- 
trated in a single pair of modes, whose amplitude is ac- 
curately described by a C-number as in the mean-field 



theory. Fluctuations give only a small correction. In par- 
ticular, the signal and idler phases are constant through- 
out the whole two-dimensional sample. As it is usual 
in the spontaneous symmetry breaking phenomena, this 
does not violate the U(l) symmetry of the system Hamil- 
tonian, since the phase is chosen at random. 

These behaviours are strictly analogous to what happens 
in a Bose gas at thermal equilibrium when the temper- 
ature is lowered across the Bose-Einstein condensation 
temperature Tbec- For T > Tbec-, the gas is a ther- 
mal, incoherent, one and is spread among many low- 
lying modes. Its coherence extends for a finite distance, 
generally smaller than the interparticle spacing, and the 
same coherence length determines the spatial behaviour 
of both the phase (first-) and the density-density (second- 
order) correlation function of the gas. For T — > Tg EC , 
the coherence length diverges&i£. For T < Tbec a 
macroscopic population is concentrated in the single state 
of lowest energy and an off-diagonal long-range order ap- 
pears in the real-space one-body density matrix. This is 
the so-called Penrose-Onsager criterion for Bose-Einstein 
condensation^. 

Also from the experimental point of view there is a close 
relation between experiments like the ones in Ref 7 which 
measure the coherence of the light emission, and experi- 
ments which address the coherence of the atomic matter 
waves. The first-order coherence of the matter wave of 
an atomic BEC has in fact been revealed by looking at 
the interference fringes formed by two atomic wavepack- 
ets extracted at different spatial positions^, while the 
density fluctuations (quantified by </ 2 ' 3 ^(0)) have been 
measured from the interaction energy^! and the three- 
body recombination rate^: both of them are quenched 
as T is lowered below Tbec- 

Although the analogy is very close, it is important to 
stress the fundamental differences between the Bose- 
Einstein condensation phenomenon of equilibrium sta- 
tistical mechanics and the appearance of spontaneous 
coherence in the parametric emission from the present 
microcavity system. As polaritons have a finite lifetime, 
they have to be continuously injected in the system from 
an external pump, and the stationary state arises as a 
sort of dynamical equilibrium between the injection and 
the dissipation. Instead of a true thermodynamical equi- 
librium, we are therefore dealing with the stationary state 
of a driven- dissipative system^. Hence, there is no rea- 
son for the critical properties to be the same, and the 
general theorems of equilibrium statistical mechanics do 
not necessarily apply. 



VII. MODULATIONAL INSTABILITY AND 
QUANTUM IMAGES 

An alternative interpretation of the parametric instabil- 
ity can be put forward in terms of modulational instabil- 
ity, the initially spatially homogeneous intensity profile 
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becoming unstable with respect to a periodic modulation 
which breaks the translational symmetry. Behaviours of 
this kind are well-known in the framework of nonlinear 
dynamics and pattern formation 3 —, and have analogies 
with the transition from the liquid to a cristalline solid 
state. The concept of modulational instability has re- 
cently played a central role in the understanding of ex- 
periments with atomic Bose-Einstein condensates^. 

The physical origin of the intensity modulation and its re- 
lation with the parametric instability can be understood 
by looking at the intensity modulation pattern formed by 
the interference of the pump, the signal and the idler. In 
k space, the pump, the signal and the idler are localized 
around respectively k p , k s and k, = 2k p — k s and the 
polariton wavefunction can be written as: 



4 P) e lk " x 



- x + VA (s) (x) e lk = x + ^ w (x)e 



-.(*)/ 



(28) 



where ip^^(pc) are slowly varying functions as compared 
to the period £ p — 27r/|k p — k s | of the intensity modu- 
lation. As a shift of the signal and idler phases by ±0 
corresponds to a spatial translation of the intensity mod- 
ulation pattern by Ax — £ p (j)/2Tr along the direction of 
k s — k p , the spontaneous breaking of the U(l) symmetry 
corresponding to the signal and idler phases can also be 
interpreted as the spontaneous breaking of the transla- 
tional symmetry along the direction of k p — k s . As the 
position of the fringes is random, the mean intensity pro- 
file remains flat, which guarantees that the translational 
invariance of the system Hamiltonian Q is preserved. 

On the other hand, fringes appear clearly in the intensity 
correlation function of the polariton field 

G^(x) = (& c (y) *^(y + x) * c (y + x) c (y)), (29) 

which is plotted in figEJ Experimentally, (x) can be 
measured by simply isolating in the near-field the light 
emitted from a pair of points y and y + x spaced of x and 
then correlating the corresponding intensity fluctuations. 
As the intensity pattern arises from the interference of 
signal, pump and idler modes, no preliminar angular se- 
lection must now be performed and, for this reason, the 
symbol G' 2 -*(x) does not carry a bar over it. 

A detailed study of the appearance of the fringes, effect 
which in quantum optics goes under the name of quantum 
image^^, was given in Ref|23for the simple case of an op- 
tical parametric oscillator formed by a medium in a 
planar cavity. Here we want to emphasize the relation be- 
tween the phase coherence properties of the signal/idler 
emission, quantified by g^(x) defined in l|24|) and the co- 
herence length of the periodic modulation of G^-'(x). As 
a shift of the signal and idler phases by ±</> corresponds 
to a spatial translation displacement of the intensity pat- 
tern, a spatial rigidity of the phase of ijj( s < 1 ) (x) transfers 
into a corresponding rigidity of the intensity modulation 
pattern. On the other hand, when the phase of ip( s ' l \x) 
is able to freely fluctuate in space, the intensity modu- 
lation pattern is washed out. This can be observed in 
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FIG. 8: Intensity correlation function G' 2 '(x) of the in- 
cavity polariton field for values of the pump frequency hS p = 
— 0.55meV (upper panel), hS p — — 0.53meV (central panel) 
and hS p = —0.52 meV (lower panel) respectively below and 
around the mean-field threshold. The cavity and the other 
pump parameters are the same as in figQ 



the Figs |S] and |§l which show G^ 2 ^(x) for the two cases 
respectively below and above the threshold: in the first 
case, the spatial extension of the modulation pattern is 
finite and equal to the phase coherence length £ c , while in 
the second case the pattern extends over the whole sys- 
tem. A related interference effect has been used in Refl26l 
to infer the coherence length of an atomic Bose-Einstein 
condensate from a measurement of intensity correlations. 

In Ref,28, the spontaneous breaking of translational 
symmetry in a planar optical parametric oscillator was 
pointed out as a slowing down of the diffusive motion 
of the modulation pattern during the time evolution of 
a single realization. Here we show another interesting 
consequence of the same phenomenon in the intensity 
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FIG. 9: Intensity correlation function G' 2 ' (x) of the in-cavity 
polariton field for a pump frequency h8 p = —0.5 meV above 
the threshold. The cavity and the other pump parameters are 
the same as in fig0 



profile of a single realization of the Wigner-Monte Carlo, 
that is in an instantanous measurement of the local in- 
tensity. While below threshold the modulation pattern is 
too weak to be detected on a single realisation and has to 
be obtained by taking the average over many realisations, 
it appears clearly although somewhat noisy in a picture 
taken above threshold when the occupation of the signal 
and idler modes is large and coherent ffig llOfl . Obviously, 
the position of the pattern is random and varies from one 
realisation to the next. 



VIII. ONE-DIMENSIONAL SYSTEM: LOCAL 
VS. GLOBAL COHERENCE 

All the physical discussions presented in the previous 
sections were based on a mean-field picture: below 
threshold, many modes are incoherently populated, while 
above threshold, there is a macroscopic occupation of 
a single signal/idler pair of mode and fluctuations give 
only a small correction. This interpretation was vali- 
dated by the numerical results of the complete Wigner- 
Monte Carlo simulations which were performed for two- 
dimensional systems of side L x = L y = 25 /im: above 
the threshold, coherence of the signal emission extended 
over the whole two-dimensional sample giving a sort of 
non-equilibrium long-range order. 

It is a well-known fact of equilibrium statistical mechan- 
ics that the spontaneous breaking of a U(l) symmetry 
and the appearance of a long-range order is possible only 
in dimension d > 2 at finite temperature and in d > 2 at 
zero temperature (Mermin- Wagner theoremiSi^) , while 
otherwise coherence extends over a finite distance only. 
In order to investigate the possibility of observing such 
behaviours in a non-equilibrium context, we have per- 
formed Wigner-Monte Carlo simulations for a system 
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FIG. 10: Intensity profile |^c(x)| 2 of a single realization 
for pump frequencies hS p = — 0.53meV (upper panel) and 
hS p — — 0.5meV (lower panel) respectively below and above 
the mean-field threshold. The cavity and the other pump 
parameters are the same as in fig0 exception made for the 
coupling constant g and the driving intensity I p respectively 
increased and decreased by a factor 5 so as to clarify the 
picture without affecting the mean-field physics. 
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FIG. 11: First- (upper panel) and second-order (lower panel) 
coherence functions <^ 1,2 '(x) of the selected signal emission 
for pump frequencies respectively h5 p = —0.53 meV below 
the mean-field threshold (solid) and hS p = —0.51 meV above 
the mean- field threshold (dotted). The cavity and the other 
pump parameters are the same as in fig0 exception made 
for the geometry: the system is here one-dimensional and 
L x = 400 /um. 
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with the same pump parameters as in the previous sec- 
tions, but in a completely different cavity geometry. In- 
stead of a square of size L x — L y — 25 /im, we consider 
here a one-dimensional system of length L x = 400 //m. 
Along this axial direction, the polaritons are free and 
periodic boundary conditions are assumed, while the 
transverse motion is assumed to be frozen in the low- 
est state by a strong lateral confinement. This kind of 
reduced-dimensionality systems are well within the pos- 
sibilities of actual semiconductor technology: effectively 
one-dimensional systems has been recently demonstrated 
by several groups^. 

The results of numerical calculations for the first- and 
second-order coherence of the signal emission are shown 
in figEJ below the threshold, the behaviour is similar 
to the one already observed in the two-dimensional case 
(figs and 0J: the phase coherence extends over a finite 
distance and intensity fluctuations due to the Hanbury- 
Brown and Twiss bunching effect are important ~ 2. 

Above the threshold, the intensity fluctuations are 
strongly suppressed g^ ~ 1 as in the two-dimensional 
case of fig[7| while the phase coherence of the signal ex- 
tends over a longer but still finite distance. The large 
distance limit of is in fact and no long-range co- 
herence exists. Similar results were discussed ir>2& for the 
case of a roll-pattern instability in an optical cavity. 

In the language of equilibrium statistical mechanics, phe- 
nomena of this kind are well-known and usually go under 
the name of quasi- condensate s 24 ^^ , so to stress the fact 
that a local order is present (signalled by the reduced 
value of the intensity fluctuations), but not the global 
one which characterizes a true Bose-Einstein condensate. 
The long range order is in fact destroyed by the long 
wavelength fluctuations of the order parameter, in our 
case the local phase of the signal emission^!. 

A final remark: recently, the realization of a BEC of 
polaritons was claimed in RefllOl starting from the ob- 
servation of a reduced value of g( 2 >. Unless a detailed 
analysis of the geometry of the system is performed, this 
observation does not appear as being a conclusive proof 
of polariton BEC, as the mere reduction of intensity fluc- 
tuations does not univocally correspond to the presence 
of a long-range order. In other terms, the possibility 
of having a quasi-condensate of polaritons has not been 
ruled out by the available experimental data. 



IX. CONCLUSIONS AND PERSPECTIVES 

In the present paper we have presented a comprehensive 
study of the coherence properties of the parametric emis- 
sion from a planar semiconductor microcavity for pump 
parameters spanning across the threshold for the para- 
metric oscillation. 

The numerical calculations have been performed by 
means of a Monte Carlo technique based on the stochastic 



differential equations of the Wigner representation of the 
coupled exciton and cavity-photon quantum fields. In- 
spired by the first experimental results reported in Ref0, 
we have calculated the coherence properties of the in- 
cavity polaritonic field and, more specifically, of the light 
which is then emitted in the signal beam. In particu- 
lar, we have analyzed the behaviour of the first and the 
second-order coherence functions when the pump param- 
eters cross the threshold for the parametric oscillation. 
The numerical results are physically understood in terms 
of a non-equilibrium phase transition occurring at the 
parametric threshold, a phase transition which closely 
ressembles to a non-equilibrium Bose-Einstein condensa- 
tion. 

Below threshold, the emission is an incoherent one, and 
the first- and second-order correlation functions show a 
finite coherence length for both the phase correlation and 
the intensity fluctuations. As the pump approaches the 
threshold, the coherence length becomes larger and larger 
and eventually becomes macroscopic at the threshold. In 
momentum space, this corresponds to a broad momen- 
tum distribution which gets narrowed into a delta-like 
peak at the threshold. 

Above threshold, the emission is instead a coherent one, 
both the signal and the idler intensities becoming of 
the same order as the pump one, and the correspond- 
ing momentum distribution being sharply peaked. For 
two-dimensional systems of size comparable to the one of 
actual experiments with planar microcavities, first- and 
second-order coherence extends in the whole system. For 
systems of reduced dimensionality such as photon wires, 
the numerical results suggest that a local order persists in 
the sense that intensity fluctuations are strongly reduced, 
but no true long-range order is present, as the global 
phase coherence is destroyed by long-wavelength phase 
fluctuations. This phenomenon is a non-equilibrium ana- 
log of the so-called quasi-condensation phenomenon of 
reduced-dimensionality Bose gases at equilibrium. 

Although the results of the present paper have been ob- 
tained for a very specific system of actual experimental 
interest, we expect the main properties of the critical 
point at the threshold to be very general for a whole 
class of non-equilibrium systems. Also the absence of 
long-range order in low-dimensional systems is expected 
to be a general fact in non-equilibrium systems: inspired 
by the Mermin- Wagner theorem of equilibrium statis- 
tical mechanics, we are actually trying to characterize 
the possibility of a long-range order in non-equilibrium 
systems of reduced dimensionality, and to determine the 
corresponding phase coherence length. In order to have 
a closer comparison with experiments, the effect of a 
disorder (which can be straightforwardly taken into ac- 
count via the potentials Vx,c( x ) m Eq^l on the non- 
equilibrium critical properties will be adressed as well. 
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The analysis of the reflected light requires more attention: 



the field emitted by the cavity through the front mirror 
has in fact to be added to the direct reflection of the pump 
beam on the front mirror. 

This can be numerically obtained by means of the same 
Wigner- Monte Carlo code by taking the mean-field limiii— 
g —* 0, E p — ► oo at a fixed g \E P \ 2 (i.e. at a fixed interaction 
energy g \ipx\ 2 )- 

In principle, a similar quasi-condensation behaviour may 
show up in two-dimensional systems as well, but the dis- 
tinction between a condensate and a quasi-condensate may 
require very exponentially large samples to rule out the 
possibility of a finite-size-induced BEG-. Analytical work 
in this direction is under way. 



